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ABSTRACT 

Many observational studies have revealed the presence of multiple stellar generations in Galactic 
globular clusters. These studies suggest that second-generation stars make up a significant fraction 
of the current mass of globular clusters, with the second-generation mass fraction ranging from ~50 
to 80% in individual clusters. In this Letter we carry out hydrodynamical simulations to explore the 
dependence of the mass of second-generation stars on the initial mass and structural parameters and 
stellar initial mass function (IMF) of the parent cluster. We then use the results of these simulations 
to estimate the fraction /sg,h of the mass of the Galactic stellar halo composed of second-generation 
stars that originated in globular clusters. We study the dependence of /sg.h on the parameters of the 
initial mass function of the Galactic globular cluster system. For a broad range of initial conditions, 
we find that the fraction of mass of the Galactic stellar halo in second-generation stars is always small, 
/SG,H < 4 - 6% for a Kroupa-1993 IMF and /sg.h < 7 - 9% for a Kroupa-2001 IMF. 
Subject headings: globular clusters: general 



1. INTRODUCTION 

During the last decade, spectroscopic and photometric 
studies have provided strong evidence for multiple stel- 
lar populations in a number of Galactic globular clusters. 
In contrast to the standard formation scenario, accord- 
ing to which globular clusters are 'simple stellar pop- 
ulations' composed of stars all with the same age and 
chemical composition, these studies have revealed star- 
to-star variations in the abundances of light elements in 
main sequence and red giant stars, suggesting a second 
stellar generation (hereafter SG) that must have formed 
out of gas processed through a first generation (hereafter 
FG) of stars. SG stars are identified by their anomalous 
abundances of Na, O, Mg and Al relative to normal halo 
field stars of similar metallicity (see e.g. Gratton et al. 
2004, Carretta et al. 2009a,b). 

Accurate photometric studies of a few globular clusters 
have strengthened the evidence for multiple populations 
by showing their fingerprints (multiple main sequences 
and sub-giant branches) in the color-magnitude diagram 
(see e.g. Lee et al. 1999, Pancino et al. 2000, Bedin et al. 
2004, Piotto et al. 2007, Milone et al. 2008, Marino et 
al. 2008, Anderson et al. 2009). The most striking main 
sequence anomalies are populations with helium abun- 
dances much larger than the standard Big-Bang value 
expected to be universally present in the oldest stars. 
The helium overabundance is interpreted as an evolu- 
tionary by-product of the progenitors of the anomalous 
SG stars. The leading models for the source gas from 
which SG stars formed involve rapidly rotating massive 
stars (e.g. Decressin et al. 2007; see also de Mink et 
al. 2009 for a more recent work suggesting massive bina- 



ries as a possible polluting source) , or intermediate-mass 
asymptotic giant branch (AGB) stars (e.g. Ventura et al. 
2001; see also Renzini 2008 for a review and references 
therein). 

Prior to the discovery of the multiple main sequences, 
the presence of stars with helium anomalies had been 
predicted based on the very extended horizontal branches 
of some clusters (D'Antona et al. 2002, D'Antona and 
Caloi 2004). Analysis of horizontal branch morpholo- 
gies of several clusters also led to the first suggestion 
that the second stellar generation could account for a 
high fraction — from 30% to 100% — of the total number 
of stars (D'Antona & Caloi 2008). While the high-quality 
photometric data necessary to detect multiple popula- 
tions have been so far obtained for only a few clusters, 
a recent spectroscopic survey of about 2000 stars in 19 
Galactic globular clusters has shown evidence of multiple 
populations in all the clusters observed, and in all cases 
SG stars account for a significant fraction (50 — 80%) of 
the cluster mass (Carretta et al. 2009a). 

D'Ercole et al. (2008) studied the formation and dy- 
namical evolution of multiple populations in globular 
clusters, and found that two distinct evolutionary pro- 
cesses can cause stars to escape and populate the Galac- 
tic stellar halo during cluster dynamical evolution. At 
early times, stars escape over the tidal boundary as the 
cluster expands due to the loss of primordial gas and FG 
SN ejecta. Since SG stars form preferentially in the cen- 
tral regions of the cluster, almost all of the escapers at 
this stage are FG stars. The result is that the SG fraction 
among bound cluster stars increases dramatically. As the 
cluster evolution continues into the phase dominated by 
two-body relaxation, the two populations mix and the 
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fractional escape rates of the FG and SG stars due to 
evaporation tend to equalize, stabilizing the SG fraction 
(see also Decressin et al. 2008 for another study of the 
mixing process). Thus the vast majority of escapers from 
the cluster are FG stars. Two recent spectroscopic stud- 
ies (Carretta et al. 2010, Marteh & Grebel 2010) have 
found that the vast majority of halo stars studied have 
abundances typical of FG stars in clusters; only about 
1.5-2.5% of the stars are Na-rich (Carretta et al. 2010) 
and CN-strong (Martell & Grebel 2010), and hence clas- 
sifiable as SG stars. 

Many general studies have addressed the evolution of 
the Galactic globular cluster system and the dispersal of 
cluster stars in the halo (see e.g. Kroupa & Boily 2002, 
Baumgardt et al. 2008, Portegies Zwart et al. 2010, Ves- 
perini 2010 and references therein). However, this Letter 
is specifically focused on SG stars. Our goals here are 
to estimate both the number of SG stars that may have 
formed in the Galactic globular cluster system and the 
fraction, /sg,h, of SG stars formed in globular clusters 
that now populate the Galactic stellar halo, and to il- 
lustrate the link between SG formation, globular cluster 
evolution and the formation history of the Galactic halo. 
In fj2l we present the results of hydrodynamical simula- 
tions exploring the dependence of the mass of SG stars 
formed on cluster initial properties. In [J3]we introduce 
our numerical framework to estimate the fraction of SG 
stars in the halo, and discuss its theoretical and obser- 
vational ingredients. In S]4] we present and discuss our 
results. 

2. MASS OF SECOND-GENERATION STARS AS A 
FUNCTION OF A CLUSTER INITIAL PROPERTIES 

In order to calculate the total mass of SG stars formed 
in a cluster as a function of cluster initial properties, 
we have carried out a series of one-dimensional hydrody- 
namical simulations modeling the SG formation process. 
The technical details of our simulations are described in 
D'Ercole et al. (2008). We outline here only the aspects 
relevant to this Letter. 

We assume that the FG cluster initially follows a King 
(1962) density profile 
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out to some truncation radius R. Stellar masses are dis- 
tributed in the range 0.1 < m/niQ < 100 according to 
either a Kroupa et al. (1993; hereafter K93) or a Kroupa 
(2001; hereafter KOI) stellar IMF. We model the forma- 
tion of SG stars from the ejecta of AGB stars and assume 
that the SG formation phase extends from about 30 Myr 
after the FG cluster forms until ^100 Myr (should the 
SG star formation epoch extend to ^ 150 Myr, and hence 
include ejecta of stars down to 4 Mq, the SG mass can be 
obtained by multiplying the values presented below by a 
factor ^ 1.4). As discussed in D'Ercole et al. (2008) and 
D'Ercole et al. (2010) in order to explain the chemical 
properties of SG stars it is necessary to invoke the pres- 
ence of pristine gas mixed in with the matter from which 
SG stars formed (see also Pfiamm-Altenburg & Kroupa 
2009 for a study of the possible accretion of gas from 
the ISM during the SG formation process). As shown by 
D'Ercole et al. (2008, 2010), a mass of pristine gas equal 



to about 50% of the total mass of SG star-forming gas 
is required to reproduce the abundance patterns (Na-0, 
Mg-Al anticorrelations and the helium distribution func- 
tion) observed in the massive cluster NGC2808 and in 
the low-mass cluster M4 (see also Ventura et al. 2009, 
Ventura & D'Antona 2009). 

However, additional systematic studies of the chemical 
properties of the individual clusters for which spectro- 
scopic data are available will be required to explore the 
allowed mass range of pristine gas and any possible de- 
pendences on protocluster conditions. Here we simply 
assume that pristine material makes up 50% of the total 
mass of gas forming the SG stars. Results for different 
pristine gas fractions can be easily obtained from those 
presented in this Letter. 

The two panels of Fig. [1] show the fractional mass of 
SG stars formed, £, = {MsG/M)init, as a function of the 
total cluster mass (M), for models with different combi- 
nations of values of the half- mass radius {Rh — 1, 2, 4 
pc), of the truncation radius {R = 20, 40, and 80 pc; 
these values encompass the range of R for the majority 
of Galactic globular clusters, for which the modal value 
is i? ~ 40 pc; see e.g. Mackey & van den Bcrgh 2005) 
and for a K93 or KOI IMF. As the cluster mass increases 
for given values of R and Rh, so too do the cluster es- 
cape speed, the fraction of AGB ejecta retained, and the 
total mass of SG stars formed. For the structural pa- 
rameters explored, only clusters initially more massive 
than ^ 10*'^ — 10^ M© can retain enough AGB ejecta to 
form an appreciable fraction of SG stars. The difference 
between the values of ^ obtained for the K93 and the 
KOI IMF is a consequence of the fact that the KOI IMF 
contains a larger fraction of the cluster mass in the mass 
range of the polluting AGB stars. 

3. SECOND-GENERATION STARS IN THE GALACTIC 
STELLAR HALO 

A complete model for the determination of the frac- 
tion of SG stars in the Galactic halo would require ac- 
curate knowledge of the initial global properties of the 
Galactic globular cluster system as well as its dynami- 
cal history, the initial internal structure and subsequent 
dynamical evolution of clusters with multiple stellar pop- 
ulations, and a comprehensive theory of how the stellar 
halo formed. Much work remains to be done on all of 
these problems. Considering the uncertainties involved 
in all of these key ingredients, we elect to proceed here 
in a simplified way, and our estimation of /sg,h proceeds 
as follows. 

For a given value of the total mass of the Galactic 
globular cluster system, Af^*^, and a given initial glob- 
ular cluster system mass function (hereafter IGCMF), 
we calculate the total mass of SG stars formed in the 
cluster system, Mg'cGCS^ using the results of the hydro- 
dynamical simulations described in the previous section. 
M'^SG^GCS m^-y be written, in general, as an integral over 
cluster structural properties 
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where f{M,Rh,R) is the joint distribution function of 
initial cluster mass and structural parameters. 
We assume either a Schechter IGCMF (see e.g. Burkert 
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Fig. 1. — Upper panel: Ratio § of the total initial mass in SG 
stars to the total initial mass of the cluster M, as a function of M, 
for models with fl = 40 pc and Rh = 1 pc, K93 IMF (filled circles), 
Rh = 2 pc, K93 IMF (asterisks), Rh = ^ pc, K93 IMF (crosses), 
i?^ = 2 pc, KOI IMF (filled triangles). Lower panel: ^ as a function 
of M , for models with Rf^ = 4 pc, K93 IMF and different values of 
the cluster truncation radius R: R = 20 pc (filled circles), R = 40 
pc (crosses), i? = 80 pc ( open circles). 
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^^^ oc M-^e-^'/^'^ for Mio^ < M < M^p, (3) 

or a simple power-law (e.g. Whitmore 2003) 

dN 

— oc M-" for Mior, < M < M^p. (4) 

With these assumptions, and for given values of Rh and 
R at the time of SG formation (30 < i < 100 Myr), Eq.[2] 
becomes 

aM,Rh,R)M-^dM. (5) 
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Fig. 2. — Ratio, rj, of the total initial mass of SG stars in the 
Galactic globular cluster system to the total initial mass of the 
cluster system as a function of the index a of the IGCMF. Differ- 
ent lines correspond to models with the following set of properties 
(IMF, Rh, R, \og(Mc) (or PL=power-law), log(Af;o„)): (upper 
panel) (KOI, 2, 40, PL, 4) upper black solid line; (KOI, 2, 40, PL, 
3) lower black solid line; (KOI, 2, 40, 6.75, 4) upper green dotted 
line; (K93, 1, 40, PL, 4) upper red dashed line; (K93, 4, 40, PL, 4) 
lower red dashed line; (K93, 2, 40, PL, 4) upper blue dot-dashed 
line; (K93, 2, 40, PL, 3) lower blue dot-dashed line; (K93, 2, 40, 
6.75, 4) lower green dotted line. (Lower panel) (K93, 4, 20, PL, 4) 
upper black solid line; (K93, 4, 20, PL, 3) lower black solid line; 
(K93, 4, 20, 6.75, 4) green dotted fine; (K93, 4, 80, PL, 4) red 
dashed line. 



Neglecting mass loss due to stellar evolution, the total 
mass of the stellar halo in the form of SG stars, either 
from dissolved globular clusters or evaporated from those 
remaining, is the difference between A^5?g*gcs' ^^'^ ^^^ 
current total mass of SG stars still in clusters, Mg^Q^g. 
A systematic observational study of SG stars in Galac- 
tic globular clusters is still lacking, and will require sig- 
nificant effort. Currently the only study affording an 
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estimate of the fraction of SG stars in a significant num- 
ber of Galactic clusters, as well as exploration of possible 
trends of the SG fraction with other cluster properties, 
is the spectroscopic survey of ^ 2000 stars in 19 Galactic 
globular clusters carried out by Carretta et al. (2009a, b). 
They find that all clusters studied host SG stars, and 
that the fraction of individual cluster masses in SG stars 
ranges from about 50% to more than 80%. Hence we 
define the current cluster SG mass fraction i^sG^cs ^y 

(6) 
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where Afg^ ~ 2 x 10' Mq is the current total mass 
of halo Galactic globular clusters (assuming M/L = 2; 
Mackey & van den Bergh 2005). 

Most of the clusters studied are more massive than 
~ IO^Mq, and it is not known whether such a large SG 
fraction is common to all Galactic globular clusters or 
only to those currently more massive than some thresh- 
old. For purposes of this Letter we simply assume as 
a reference value -F'sg'gcs — O-^- Results for different 
values of -Fsg^cs '^^^ ^^ very easily derived from those 
shown in this Letter. 

To complete our calculation we need an estimate of 
^GCS' ^^^'^ initial total mass of all stars in globular clus- 
ters. This is not easily determined from observations of 
the current Galactic globular cluster system. Here again, 
we simply parametrize the normalization of Eq.[3]or|4]by 
defining the parameter $ = MQ-^g/Mhaio where we take 
Mfiaio = IO^Mq (see e.g. Freeman & Bland-Hawthorn 
2002). 

The current fraction of SG stars in the stellar halo, 
/sG,H: then is 



/sG,H 



]\ jinit 



]\/fnow 
^^^SG,GCS 



M, 



halo 



(7) 



Finally, if we define rj as the ratio of the total initial mass 
of SG stars in the Galactic globular cluster system to the 
total initial mass of the cluster system 

^ ^ M^acs/M}^s, (8) 

and use the definition of ^sg'gcs ^^^ Eq.[6] along with 
the values of Mhaio and Mq^c; adopted above, we can 
rewrite Eq.[7]as 



/sG,H = ^*-0.02Fg"<S- s- 



(9) 



The form of Eq.[5] allows to easily identify the different 
theoretical and observational ingredients needed for the 
determination of /sg.h- We summarize them here. 

• rj, the ratio of the total initial mass of SG stars 
to the total initial mass of the Galactic globu- 
lar cluster system, incorporates the dependency 
on the SG formation model (through the function 
^{M,Rh,R)), as well as on the stellar IMF, the 
cluster mass and structural parameter distribution 
functions (see Eq.[5] above). 

• <&, the ratio of the total initial mass of the Galactic 
globular cluster system to the current halo mass, is 
not known. The simple linear scaling of /sg,h with 
this parameter allows to easily explore the depen- 
dence of /sg.h on the total initial mass of the 
Galactic globular cluster system. 




Fig. 3. — /sG.H versus <I>. A power-law IGCMF is assumed with 
a = 1.8. Different lines correspond to the following set of properties 
(IMF, Rh, R, \og{Mia^)): (KOI, 2, 40, 4) upper black solid line; 
(KOI, 2, 40, 3) lower black solid line; (K93, 4, 20, 4) upper blue 
dashed line; (K93, 4, 80, 4) lower blue dashed line; (K93, 1, 40, 4) 
upper red dot-dashed line; (K93, 2, 40, 4) middle red dot-dashed 
line; (K93, 4, 40, 4) lower red dot-dashed line; (K93, 2, 40, 3) green 
dotted line (almost overlapping the lower blue dashed line). 



• ^SG.GCS' 



the fraction of the total current mass 
of the halo cluster population composed of SG 
stars, is a parameter to be determined observa- 
tionally. In our study, based on the spectroscopic 
survey of Carretta et aZ. (2009a, b), we adopt as a 
reference value ^sg'gcs ~ '^■^- However, in this 
case too, the simple linear scaling of /sg,h with 
^sg'gcs allows to easily calculate /sg,h for other 
values of Fgo^^cs • 

4. RESULTS AND DISCUSSION 

In Fig. [5] we show the dependence of rj on the power- 
law index a of the IGCMF. Different lines correspond to 
different pairs of values of Rt and i?, a K93 or a KOI IMF 
and different values for the parameters of the IGCMF 
{Mc and Af/o^; Mup is fixed at lO'' Mq). 

Comparison of the results for models with different val- 
ues of R shows no significant differences. 

An increase in Mc and/or Miow increases the fraction 
of the total initial mass in clusters massive enough to 
form SG stars (M > 10"-^ - 10^; see Fig. [J), and hence 
increases rj. Similarly, as a decreases the fraction of total 
mass in massive clusters increases, and so does rj. 

For a fixed value of i?, decreasing the value of Rh 
adopted has the obvious consequence of increasing the 
fraction of AGB ejecta retained and the amount of SG 
stars formed in lower-mass clusters. The absolute upper 
hmit {rj^ax ^ 0.08 for a K93 IMF and rj^^ax ^ 0.125 
for a KOI IMF) would be attained by assuming (un- 
realistically) that all clusters (including the low-mass 
ones) were compact enough at the time of SG forma- 
tion (30 < t < 100 Myr) to retain aU the AGB ejecta and 
reach a fraction of SG stars formed ^ ~ 0.08 or ^ ~ 0.125 
for a K93 and a KOI IMF respectively. 
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From the results shown in Figj5] one can easily calcu- 
late /sG,H as a function of $ and -FgQ'QQg for any of the 
combinations of IGCMF, cluster structural parameters 
and IMF explored. For example, Fig.[3l shows the de- 
pendence of /sG,H on $ for a power-law IGCMF with 
a = 1.8, for various values of Mio^^ R, Rh and for ei- 
ther a K93 or a KOI IMF (and, as discussed above, we 
have adopted -F'sg'gcs ~ ^.5). If clusters were at the 
time of SG formation (30 ^ t ^ 100 Myr) more compact 
than assumed here, the lines of /sg,h vs $ would fall 
between those shown in Fig. [3] and that corresponding to 
the absolute upper limit easily obtained using Eql9] and 
the maximum value of rj {rimax — 0.08 for a K93 IMF 
and n^ax ^ 0.125 for a KOI IMF). 

Figs. [2] and [3] demonstrate that, for the broad range 
of values of the IGCMF parameters considered, the pre- 
dicted fraction of SG stars in the halo is always small, 
/SG,H < 4 - 6% for a K93 IMF and /sg,h < 7-9% 
for a KOI IMF . This result is consistent with the small 
fraction (■~1.5-2.5%) of halo stars identified as SG stars 
by Carretta et al. (2010, based on the stars' sodium 
abundance) and Martell and Grebel (2010, based on GN 
and CH bandstrengths). For comparison, 50% of glob- 
ular cluster stars are identified as SG when the same 



classification criteria are adopted (Kraft 1994, Carretta 
et al. 2010, MarteU & Grebel 2010). We note that for 
these estimates of /sg.Hi using, for example, a power- law 
IGCMF with a — 1.8 as a reference model, our calcula- 
tions (Fig.El) with a KOI IMF (K93 IMF) imply that a 
fraction ranging from about 20% to about 40% (30% to 
about 60%) of the Galactic stellar halo must be composed 
of stars originally formed in globular clusters. Although 
this estimate is very uncertain and both more calcula- 
tions and further observational studies of SG stars in 
the halo will be required to refine it, it is interesting to 
note how the fractions of SG stars in globular clusters 
and in the Galactic stellar halo connect and constrain 
both models for the formation and dynamical evolution 
of multiple stellar generations in globular clusters and 
the formation history of the Galactic halo. 
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